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'  ASSTRACT 

We  formulate  a  mathematical  model  for  the  mean  disturbance  created  in  a  linearly  strat¬ 
ified  fluid  by  the  growth  and  collapse  of  a  localized  region  of  turbulent,  mixed  fluid.  Linear¬ 
ized  differential  equations  for  the  disturbance  are  obtained  from  the  Navier-Stokes  equations 
under  the  Boussinesq  approximation  after  various  intuitive  and  empirical  considerations.  The 
divergence  of  the  mean  turbulent  density  flux  enters  these  equations  as  the  sole  forcing 
function  for  the  disturbance*  this  differs  from  previous  theoretical  models  in  which  an  initial, 
perturbed  density  field  acts  as  the  forcing  function. 

We  replace  'the  original  expression  we  find  for  the  forcing  function  by  an  approxima¬ 
tion  devised  to  simplify  th£  numerical  evaluation  of  the  disturbance.  We  then  obtain  a 
solution  for  the  disturbance  in  a  fluid  of  constant  depth  and  compare  the  theoretical 
rwnlt»  with  measured  data.  The  comparison  indicates  that  the  theory  provides  satisfactory 
results  at  large  distances  from  the  mixed  region  but  has  deficiencies  near  the  mixed 
region.  Specific  proposals  are  made  for  correcting  the  deficiencies  of  the  model  at  small 
distances  from  the  mixed  region. 
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ABSTRACT 


We  formulate  a  mathematical  model  for  the  mean  disturbance 
created  in  a  linearly  stratified  fluid  by  the  growth  and  collapse  of  a 
localized  region  of  turbulent,  mixed  fluid.  Linearized  differential 
equations  for  the  disturbance  are  obtained  from  the  Navier-Stokes 
equations  under  the  Boussinesq  approximation  after  various  intuitive 
and  empirical  considerations.  The  divergence  of  the  mean  turbulent 
density  flux  enters  these  equations  as  the  sole  forcing  function  for 
the  disturbance;  this  differs  from  previous  theoretical  models  in 
which  an  initial,  perturbed  density  field  acts  as  the  forcing  function. 

We  replace  the  original  expression  we  find  for  the  forcing  func¬ 
tion  by  an  approximation  devised  to  simplify  the  numerical  evalua¬ 
tion  of  the  disturbance.  We  then  obtain  a  solution  for  the  disturbance 
in  a  fluid  of  constant  depth  and  compare  the  theoretical  results  with 
measured  data.  The  comparison  indicates  that  the  theory  provides 
satisfactory  results  at  large  distances  from  the  mixed  region  but  has 
deficiencies  near  the  mixed  region.  Specific  proposals  are  made  for 
correcting  the  deficiencies  of  the  model  at  small  distances  from  the 
mixed  region. 
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auxiliary  quantity 
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pressure 

fluctuating  component  of  pressure 
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pressure  in  undisturbed  fluid 
radius  of  mixed  region 
effective  radius  of  mixer 
Reynolds  stresses,  -po<u'u') 
viscous  stress  tensor 
ensemble  average  temperature 
temperature  prior  to  mixing 
time 

time  at  which  mixed  region  has  radius  r 
time  for  mixed  region  to  reach  point  (x,y) 

Brunt-Vaisala  period 

excitation  time;  interval  during  which  mixer  agitates  the  fluid 

duration  of  mixing  phase 

reference  time,  virtual  origin  of  mixed  region 

fluid  velocity 

fluctuating  component  of  velocity 
ensemble  average  velocity 
vertical  mean  velocity  component 
Fourier  transform  of  w 

vertical  velocity  component  associated  with  frequencies  less  than  N 
vertical  velocity  component  associated  with  frequencies  greater  than  N 
horizontal  coordinate 

vertical  coordinate  referred  to  center  of  mixed  region,  z  -  z0 
vertical  coordinate 


iii 


z0  depth  to  center  of  mixer 

orn  auxiliary  quantity 
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0„  auxiliary  quantity 

AH  average  temperature  duration 

AT  ensemble  average  temperature  change 
5(x)  Dirac  delta  function 

6 p  density  change  produced  by  mixing 
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p  density 
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<t>n  set  of  orthogonal  functions 

cj  frequency,  rad/s 


iv 


INTERNAL  WAVE  GENERATION  CAUSED  BY  THE  GROWTH 
AND  COLLAPSE  OF  A  MIXED  REGION 


INTRODUCTION 

The  phenomenon  of  internal  wave  generation  by  the  turbulent  wake  of  a  self- 
propelled  body  moving  through  a  density-stratified  fluid  was  first  demonstrated  by 
Schooley  and  Stewart  (1)  in  1963.  Their  procedure  involved  an  experimental  study  of 
the  disturbance  produced  by  a  small-scale,  self-propelled  body  in  a  laboratory  tank  under 
carefully  controlled  conditions. 

Schooley  and  Stewart  noted  that  the  horizontal  and  vertical  dimensions  of  the  tur¬ 
bulent  wake  initially  increased,  but  eventually  the  vertical  dimension  reached  a  maximum 
and  afterward  decreased  in  value.  They  attributed  this  collapse  of  the  wake  to  the  return 
of  mixed  fluid  in  the  wake  to  an  appropriate  equilibrium  level  under  the  action  of  grav¬ 
ity.  They  also  argued  that  for  bodies  moving  at  moderate  speeds  the  collapse  process 
should  be  approximately  a  two-dimensional  process.  This  concept  of  collapse  and  the 
resulting  internal  wave  generation  as  essentially  a  two-dimensional  process  led  to  further 
experimental  work  (2-6)  in  which  the  mechanism  of  inducing  turbulence  was  by  means  of 
various  cylindrical  mixers,  chosen  so  as  to  produce  motions  whose  average  characteristics 
could  be  regarded  as  two  dimensional. 

These  later  experiments  fall  into  two  distinct  categories.  The  first  category  is  repre¬ 
sented  by  Wu’s  experiment  (2),  in  which  a  portion  of  the  fluid  is  kept  separate  from  the 
surrounding  fluid  by  a  rigid  partition  as  it  is  being  mixed;  once  turbulence  has  subsided 
in  the  mixed  fluid  the  rigid  partition  is  removed  and  collapse  ensues.  Mei  (7)  presented 
a  theoretical  analysis  of  the  collapse  phase,  and  Wessel  (8)  developed  a  numerical  method 
for  studying  the  process  based  on  the  Navier-Stokes  equations. 

The  second  category  provides  a  more  realistic  representation  of  the  interaction 
between  the  turbulent  wake  of  a  self-propelled  body  and  the  ambient  stratified  fluid  than 
does  the  first  category;  it  consists  of  those  experiments  in  which  the  mixed,  turbulent 
fluid  is  always  in  contact  with,  and  interacting  with,  the  surrounding  unmixed  fluid. 

Most  experimental  investigations  of  two-dimensional  collapse  fall  into  this  category  (see 
Refs.  3-6).  A  characteristic  feature  of  these  experiments,  not  present  in  Wu’s  experiment, 
is  that  the  mixed  region  initially  expands  in  extent  until  some  characteristic  time  is 
reached  when  collapse  of  the  region  takes  place. 

Recently,  Schooley  and  Hughes  (9)  developed  a  theoretical  model  for  wave  genera¬ 
tion  by  collapse  and  made  some  comparison  with  the  results  of  an  experiment  of  the 
second  category.  They  found  rough  agreement  with  measurements  made  at  large  dis¬ 
tances  from  the  mixer;  however,  they  took  their  origin  of  time  as  arbitrary  and  deter¬ 
mined  it  by  requiring  as  good  a  match  as  possible  between  their  results  and  the  measured 
data.  An  undesirable  feature  of  their  model,  which  explains  their  need  for  an  arbitrary 
time  origin,  is  the  neglect  of  the  initial  growth  stage  of  the  mixed  region  and  its  role  in 
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the  wave  generation  process.  They  assume,  in  fact,  that  wave  generation  begins  essen¬ 
tially  at  the  same  time  as  does  collapse  of  the  mixed  region.  This  assumption  is  not  in 
accord  with  the  experimental  results  of  Williams  (6)  who  made  measurements  near  the 
mixer  and  found  that  wave  generation  begins  as  soon  as  the  mixer  is  activated. 

In  this  report  we  develop  a  theoretical  model  for  internal  wave  generation  which 
takes  approximate  account  of  the  entire  history  of  the  mixed  region.  We  do  this  by 
specifying  on  empirical  and  intuitive  grounds  only  those  properties  of  the  mixed  region 
which  we  feel  are  important  for  the  process  of  wave  generation.  Several  different  models 
could  be  developed  by  this  procedure,  and  we  have  selected  one  of  the  simplest  as  a  can¬ 
didate  for  study. 

We  base  our  model  on  a  linearized  version  of  the  Navier-Stokes  equations  in  the 
form  applicable  to  the  ensemble  average  motion.  The  linearization  is  motivated  by  prac¬ 
tical  considerations  rather  than  by  experimental  evidence:  Linearizing  the  Navier-Stokes 
equations  produces  a  system  of  equations  amenable  to  mathematical  analysis  whose  appli¬ 
cability  to  the  physical  process  of  interest  can  then  be  examined. 

Our  linearized  equations  contain  the  mean  turbulent  density  flux  and  the  Reynolds 
stresses  as  forcing  terms  for  the  disturbance.  None  of  these  quantities  have  been  mea¬ 
sured  in  the  experiments  cited  above,  and  little  appears  to  be  known  about  their  behavior 
in  a  stratified  fluid.  On  intuitive  grounds,  we  take  the  mean  turbulent  density  flux  as  the 
primary  forcing  term  and  obtain  the  information  about  this  flux  required  to  complete  the 
mathematical  model  through  empirical  considerations.  We  then  solve  the  differential 
equations  of  the  model  for  the  particular  case  of  a  linearly  stratified  fluid  of  constant 
depth.  The  analytical  solution  is  compared  to  experimental  results  of  Schooley  and 
Hughes  (9)  which  apply  to  points  at  some  distance  from  the  mixer.  We  also  compare  the 
theory  with  measurements  made  by  Williams  (6)  near  the  mixed  region.  In  this  case, 
interpretation  of  the  comparison  is  difficult  because  expeiimental  conditions  do  not  agree 
wiin  requirements  of  the  theory.  We  do  find  that  similar  trends  occur  in  both  the  theo¬ 
retical  results  and  the  measured  data.  However,  the  main  implication  of  our  comparison 
with  the  Williams  data  seems  to  be  that  the  present  model  is  deficient  for  points  near  the 
mixed  region.  Finally,  we  suggest  modifications  of  the  model  which  may  remove  the  de¬ 
ficiencies  in  the  solution  for  points  near  the  mixer. 


THE  MATHEMATICAL  MODEL 
Approximate  Linearized  Equations 

Our  main  objective  is  the  development  ol'  a  mathematical  model  to  describe  the  in¬ 
ternal  wave  system  created  by  the  growth  and  collapse  of  a  localized  region  of  turbulence 
in  a  linearly  stratified  fluid.f  We  assume  that  the  stratified  fluid  has  a  finite  depth  D  but 
is  effectively  infinite  in  horizontal  extent.  A  cylindrical  mechanical  device  called  a  mixer 
with  an  effective  radius  r0  and  located  at  some  depth  z0  creates  turbulence  in  the  fluid 
when  oscillated  or  rotated  back-and-forth  around  its  axis  at  a  sufficient  rate.  Figure  1 
shows  the  geometry  of  the  flow.  We  assume  that  the  mixer  acts  for  a  finite  interval  of 
time  and  produces  a  flow  whose  ensemble  average  characteristics  are  two  dimensional, 
dependent  only  on  x,  z,  and  t.  Figure  1  shows  the  coordinate  system  used  to  describe 


tOur  approach  can  be  easily  extended  to  include  more  general  stratifications,  but  the  numerical  problem 
of  evaluating  the  disturbance  becomes  complicated. 
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Fig.  1— Coordinate  system  for  description  of  flow;  rQ  designates 
the  effective  radius  of  the  mixer  and  zo  the  depth  of  the  center 
of  the  mixer 


the  flow.  It  will  be  further  assumed  that  the  structural  details  of  the  mixer  can  be  ig¬ 
nored  to  a  first  approximation.  (Inclusion  of  the  mixer  in  the  mathematical  problem  as  a 
rigid  boundary  to  the  flow  generally  leads  to  an  intractable  problem.)  In  this  section  we 
formulate  approximate  equations  for  the  description  of  internal  wave  generation  caused 
by  the  mixing  of  the  fluid. 

The  basic  equations  describing  the  flow  are  the  Navier-Stokes  equations  modified  by 
the  Boussinesq  approximation  (see  Ref.  10).  These  equations  consist  of  the  incompressi¬ 
bility  equation 


for  the  density  p,  the  solenoidal 


for  the  fluid  velocity  u,  and  the 

Po 


condition 

V  •  u  =  0 
momentum  equation 

“  -  -Vp  +  V  •  s  +  pg. 


(1) 

(2) 


in  which  p0  denotes  the  mean  density  of  the  fluid,  p  the  pressure,  S  the  viscous  stress 
tenor,  and  g  the  acceleration  of  gravity.  For  those  cases  of  inte-est  to  us  the  coefficient 
of  viscosity  p  can  be  regarded  as  constant  and  this  allows  the  momentum  equation  to  be 
written  as 


du 

dt 


-  —  Vp  +  —  9  +  i>V2u, 
Po  Po 


(3) 


with  the  kinematic  viscosity  v  =  p/p0. 


For  turbulent  flows  it  is  not  possible  to  describe  mathematically  particular  realizations 
of  the  flow;  instead  we  must  use  a  statistical  description  based  on  the  average  flow  charac¬ 
teristics.  In  our  problem  we  can  define  the  average,  or  mean,  flow  as  the  ensemble 
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average  taken  over  many  individual  realizations  of  the  flow,  each  of  which  possesses  the 
same  values  for  controllable  parameters.  We  require  then  differential  equations  applicable 
to  the  mean  flow. 

Each  flow  quantity  can  be  written  as  the  sum  of  its  ensemble  average  and  a  fluctu¬ 
ating  term  whose  mean  value  vanishes.  Thus,  we  write 

P  «  <p>  +  p' 

u  =  v  +  u' 

p  *  P  +  p\  (4) 

in  which  the  prime  quantities  denote  fluctuating  terms.  We  note  that  the  mean  flow 
quantities  generally  depend  on  time.  Differential  equations  for  the  mean  flow  result  from 
applying  the  ensemble  average  to  Eqs.  (l)-(3).  Equation  (1)  provides 

~<p)  +  v  •  V<p>  =  -V  •  D,  (5) 

in  which  D  =  <p'u’>  is  the  average  density  flux  due  to  the  fluctuating  motion.  From  Eq. 
(2)  we  obtain 


V  •  v  =  0,  (6) 

which  expresses  the  fact  that  the  mean  flow  is  solenoidal.  Equation  (3)  leads  to 

~~  +  v  •  W  =  -  -  VP  +  — g  +  uV2V  +  —V  -  R  (7) 

ot  Pq  p0  Pq 

in  which  R  =  -p0<u'u'>  designates  the  Reynolds  stresses  acting  on  the  mean  flow  as  a 
consequence  ot  the  fluctuating  velocity  field.  We  note  that  measurements  of  the 
Reynolds  stresses  and  the  average  density  flux  do  not  exist  for  those  situations  of  interest 
to  us. 

To  simplify  Eqs.  (5)-(7),  we  can  use  the  fact  that  the  Reynolds  stresses  in  Eq.  (7) 
can  be  written  as  a  pressure  term  plus  a  term  containing  sheai  stresses  alone.  We  can 
combine  the  pressure  term  with  P  to  essentially  redefine  P.  As  for  the  shear  stresses,  we 
assume  they  are  unimportant  for  describing  the  wave  generation  process.  We  note  that 
for  locally  isotropic  turbulence  these  shear  stresses  vanish.  However,  in  the  present  prob¬ 
lem  they  probably  are  important  at  least  in  the  vicinity  of  the  mixer.  We  also  assume 
that  the  viscous  stresses  do  not  play  a  crucial  role  for  the  wave  generation  process  and 
drop  them  from  Eq.  (7). 

We  write 


<p>  =  p*(z)  +  fi 

P  -  P*(z)  +  p,  (8) 


in  which  p*  and  P*  refer  to  the  equilibrium  state  of  the  stratified  fluid  prior  to  mixing 
and  satisfy 


t rfW?'**  f—-  ■** 
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dP* 

dT  =  p  *'  (9) 

We  always  have  the  condition  that  p  is  small  in  comparison  with  p*(z).  However  we 
shall  assume  that  SJP  can  be  neglected  in  com  par  son  with  Vp*  and  that  the  mean  flow 
velocity  v  can  be  regarded  as  small.  These  assumptions  do  not  stem  from  experimental 
evidence.  Rather,  they  are  made  for  the  sole  purpose  of  obtaining  an  analytically  trac¬ 
table  system  of  equations  whose  applicability  to  the  physical  process  of  interest  can  then 
be  studied.  We  should  point  out  that,  even  though  these  assumptions  may  turn  out  to  be 
untenable,  the  simplified  equations  they  produce  may  still  provide  an  approximate  de¬ 
scription  of  some  features  of  the  physical  process. 

Substituting  Eqs.  (8)  and  (9)  into  Eqs.  (5 )-(7 )  and  neglecting  those  terms  which  are 
of  second  order  in  small  quantities,  we  obtain  the  following  linear  equations  for  the  mean 
disturbance: 


and 


V  •  v  =  0 


(10) 


M 

dt 


+ 


-V  •  D, 


(11) 


in  which  w  is  the  vertical  velocity  component,  and 


dv 

dt 


-  -Vp  + 
P0 


P  0 


(12) 


Equations  (10)-(12)  represent  the  basic  differential  equations  we  shall  employ  in  de¬ 
scribing  the  internal  wave  system  produced  by  the  mixer.  They  have  been  established 
mainly  on  intuitive  grounds  because  of  the  lack  of  detailed  experimental  measurements. 


An  auxiliary  quantity  of  use  in  describing  the  disturbance  is  the  displacement  tj  de¬ 
fined  by 


3i7 

aT  =  * 


and 


T)(x,z,  0)  =  0. 


(13) 


The  function  r}(x,  z,  t)  gives  the  displacement  in  the  vertical  direction  of  the  particle 
initially  (that  is,  prior  to  mixing)  at  (x,  z)  produced  by  the  mean  motion.  In  those 
regions  of  the  flow  where  mixing  or  other  diffusive  mechanisms  do  not  occur,  the  tem¬ 
perature  of  a  fluid  particle  is  conserved  throughout  its  motion:  this  property  leads  to  a 
relation  between  the  temperature  change  at  a  point  AT  *  T  -  T*  ( T *  being  the  temper¬ 
ature  prior  to  the  start  of  mixing)  and  the  displacement  77  given  by 


AT  - 


(14) 
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This  relationship  is  quite  important  for  comparing  theoretical  results  with  experimental 
measurements  because  internal  waves  are  often  measured  through  the  time-dependent 
temperature  change  they  produce  at  a  point. 

The  form  of  Eqs.  (10)-(12)  indicates  that  the  mean  density  flux  D  acts  as  the 
forcing  function  for  the  disturbance.  This  can  be  made  more  explicit  by  integrating  Eq. 
(11)  and  using  Eq.  (13)  as  well  as  the  initial  condition  that  no  disturbance  is  present  at 
f  =  0  to  get 


D)dt. 


(15) 


This  equation  states  that  the  density  variation  at  a  point  is  the  sum  of  two  effects.  The 
first  effect  is  the  mean  convection  of  the  isopycnals  represented  by  the  first  term  on  the 
right-hand  side  of  Eq.  (15).  The  second  effect  leading  to  a  change  in  density  is  the 
presence  of  turbulence  or  mixing,  which  changes  the  mean  value  of  the  density.  This 
effect  is  represented  by  the  second  term  on  the  right  side  of  Eq.  (15).  Since  the  buoy¬ 
ancy  force  is  given  by  fig  we  can  also  consider  the  buoyancy  force  as  the  sum  of  two 
components.  One  component  is  the  reaction  of  the  stratified  fluid  to  displacement  of 
the  isopycnals  by  the  mean  flow,  which  is  always  in  a  direction  to  restore  equilibrium. 
The  second  component  of  the  buoyancy  force  is 


9 


D  )dt 


(16) 


and  depends  entirely  on  the  nature  of  the  turbulence.  Immediately  after  the  start  of 
mixing  this  second  component  of  the  buoyancy  force  begins  to  act  on  the  fluid  and 
creates  a  disturbance  which  ultimately  propagates  throughout  the  entire  fluid. 

To  complete  the  statement  of  the  mathematical  problem  we  must  specify  V  •  D  as 
well  as  the  boundary  conditions  and  initial  conditions.  Measurements  of  V  *  D  do  not 
exist  for  the  situations  of  interest  to  us,  and  we  are  forced  to  rely  on  empirical  consid¬ 
erations  in  order  to  define  a  functional  form  for  V  •  D. 


Specification  of  the  Forcing  Function 

We  find  a  functional  form  for  the  quantity  V  •  D  which  describes  the  net  mixing  of 
the  fluid  by  use  of  experimental  results  and  intuitive  arguments.  A  few  properties  of  this 
function  are  obvious:  outside  the  mixed  region  it  has  the  value  zero;  once  the  generation 
of  turbulence  by  the  mixer  has  ceased  there  is  a  finite  interval  of  time  during  which  the 
turbulence  loses  energy  by  working  against  the  stratification  and,  at  the  end  of  this  time 
interval,  the  net  mixing  is  zero,  for  all  practical  purposes.  We  designate  by  te,  called  the 
excitation  time,  the  time  interval  during  which  the  mixer  generates  turbulence  and  by  tm 
the  duration  of  the  mixing  phase,  which  is  the  time  interval  from  the  moment  the  mixer 
is  activated  to  the  time  when  net  mixing  vanishes.  Turbulence  may  still  be  present  after 
the  end  of  the  mixing  phase,  but  it  will  produce  no  net  mixing  of  the  fluid. 


The  quantity  V  *  D  will  therefore  be  nonzero  inside  the  boundary  of  the  mixed 
region  over  the  time  interval  0  <  t  <  tm .  Figure  2a  shows  a  representative  history  of  the 


WKJTM  (cm) 
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horizontal  and  vertical  dimensions  of  the  mixed  region.  The  dramatic  decrease  of  the 
vertical  dimension  while  the  horizontal  dimension  continues  to  increase  is  clearly  indi¬ 
cated.  Figure  2b  gives  the  time  variation  of  the  area  of  the  mixed  region;  the  area  was 
calculated  by  assuming  an  elliptical  cross  section  for  the  mixed  region  and  using  the 
results  in  Fig.  2a.  The  important  features  of  the  behavior  are  the  rapid  initial  increase 
and  the  attainment  of  a  constant  value.  The  increase  in  area  at  later  times  is  difficult  to 
explain  but  apparently  may  be  due  to  the  action  of  boundary  layers  on  the  side  walls  of 
the  tank  (the  side  walls  were  about  1  cm  apart).  We  do  not  consider  this  increase  in  area 
at  the  latest  times  as  playing  a  crucial  role  in  the  process  of  wave  generation.  Thus,  we 
idealize  the  behavior  by  interpreting  it  as  an  initial  increase  to  a  final,  constant  value. 


(a)  Total  vertical  and  horizontal  di¬ 
mensions  of  the  region 


(b)  Area  of  turbulent  region  determined  by  use 
of  (a)  and  the  assumption  of  an  elliptical  cross 
section  for  the  region  (based  on  Schooley  (5)) 


Fig.  2— Representative  time  history  of  turbulent  region.  Conditions  are  Brunt-Vtiisala  period  lBV  -7  s, 
mixer  radius  tq  =  0.64  cm,  excitation  time  te  =  2.5  s 


A  constant  value  for  the  area  of  the  mixed  fluid  indicates  that  entrainment  of  fluid 
has  ceased;  presumably,  the  turbulent  motion  no  longer  has  sufficient  energy  to  mix  addi¬ 
tional  fluid.  We  therefore  assume  that  mixing  is  complete  when  the  area  of  the  mixed 
region  no  longer  increases  with  time.  This  offers  a  means  for  determining  tm ,  the  time 
interval  during  which  mixing  occurs.  We  take  tm  to  depend  on  the  excitation  time  te 
and  the  Brunt-Vaisala  period  according  to  the  simple  law 


tm  =  constant  x  tBV  +  te. 

From  Fig.  2b  the  va’ue  of  the  constant  is  0.2,  which  gives 

‘m  “  0.2f BV  +  te- 


(17) 


Intuitively  we  would  expect  a  relationship  of  the  form  given  by  Eq.  (17),  since  it  states 
that  mixing  occurs  during  the  excitation  interval  plus  an  interval  of  time  dependent  on 
the  Brunt- VSisila  period.  Equation  (17)  also  seems  to  be  in  accord  with  the  measure¬ 
ments  of  van  de  Watering  (3). 
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We  now  know  the  period  of  time  over  which  V  •  D  is  nonzero.  We  further  require 
knowledge  of  the  spatial  region  in  which  this  quantity  is  nonvanishing.  To  satisfy  this 
requirement  we  assume  that  a  boundary  can  be  defined  for  the  mixed  region  which  has  a 
circular  form  over  the  interval  0  <  t  <  tm  with  a  time-dependent  radius  r(t).  This  approx¬ 
imately  agrees  with  measured  results  (see  Fig.  2a).  For  analytical  convenience  we  assume 
that 


r(t)  =  Ar0j/r*~%  forO  <  t  <  tm,  (18) 

in  which  t  =  -  t0  gives  the  virtual  origin  of  the  mixed  region.  The  quantities  A  and  f0 
can  be  evaluated  by  requiring  that  the  area  i rr2(t)  should  agree  with  measured  data.  For 
those  cases  of  interest  to  us  f0  is  about  0.1  s  and  quite  small  compared  to  the  Brunt- 
Vaisala  period.  Thus,  we  shall  shift  our  time  origin  by  t0  and  regard  the  mixed  region  as 
expanding  from  a  point.  This  allows  us  to  write 

r(f)  =  Ar0f/7  forO  <  t  <  tm,  (19) 

In  Fig.  2b  it  is  shown  that  a  2/3  law  provides  an  adequate  fit  to  the  data. 


We  expect  that  the  quantity  A  appearing  in  Eqs.  (18)  and  (19)  depends  on  the  effec¬ 
tive  frequency  of  oscillation  of  the  mixer  and  on  the  geometrical  structure  of  the  mixer. 
Measured  data  do  not  permit  this  dependence  to  be  determined.  We  shall  treat  A  as  a 
constant  to  which  a  value  will  be  assigned  on  the  basis  of  measured  results  for  each  par¬ 
ticular  experiment. 


The  main  assumption  we  use  to  obtain  a  value  for  V  •  D  is  that  mixing,  in  the  ab¬ 
sence  of  mean  motion,  produces  a  constant  degree  of  mixing  0  inside  the  mixed  region. 
The  degree  of  mixing  0  is  defined  as 


0  = 


6p 

P*(zq)  -  P*(z)’ 


(20) 


in  which  5 p  designates  the  density  change  produced  by  mixing.  Since  p*(z)  has  a  linear 
variation  we  can  write  Eq.  (20)  as 


5p  =  0  ^  (z°  -  z)  (21) 

inside  the  mixed  region.  Figure  3  shows  the  density  variation  along  a  vertical  line 
through  the  center  of  the  mixed  region. 

Using  Eqs.  (15)  and  (19),  and  the  assumption  of  a  constant  degree  of  mixing  in  the 
absence  of  mean  motion,  we  find  that  the  mean  turbulent  density  flux  must  satisfy 

-v  •  D  =  0^— -(z o  "  *)$(*  -  ffr(r))W(tm  -  t)  fort  >  0,  (22) 

in  which  8(x)  designates  the  delta  function  and  H(x)  the  unit  step  function,  with  the 
properties 


j 
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;-jKWfn  y*  i»fcw^= 


DENSITY 


INITIAL  DENSITY  VARIATION 
\  WITH  SLOPE 


r(t)  t  \  ALTERED  DENSITY  IN  MIXED  # 
\  REGION  WITH  SLOPE  (1-/9)  4^ 

8/)-4w-  dI 


Fig.  3— Density  variation  along  line  through 
center  of  mixed  region;  20  gives  depth  of  mixed 
region,  r(t)  the  radius  of  mixed  region  at  time 
t,  0  the  degree  of  mixing,  and  6 p  the  density 
change  produced  by  mixing  in  the  absence  of 
mean  motion 


H(x )  =  1  for*  >  0, 
H(x)  =  0  for*  <  0. 


The  quantity  tb(r)  gives  the  time  at  which  the  boundary  of  the  mixed  region  reaches  a 
point  at  a  distance  r  from  the  center  of  the  mixed  region;  analytically. 


(*)' 


We  now  have  an  explicit  functional  expression  for  the  forcing  function.  However, 
when  the  differential  equations  are  solved  using  this  expression,  we  obtain  a  solution  for 
the  disturbance  in  which  four  integrations  appear.  This  leads  to  a  complicated  numerical 
problem  which  we  prefer  to  avoid.  Thus,  instead  of  employing  the  form  of  V  •  D  given 
by  Eq.  (22),  we  shall  use  an  alternative  expression  for  V  •  D  which  approximately  pro¬ 
duces  the  same  effects  as  the  expression  in  Eq.  (22).  We  take  this  approach  for  the  sake 
of  mathematical  convenience  and  note  that  several  approximate  forms  for  V  •  D  could  be 
developed.  It  is  difficult  to  select  the  best  approximate  form  to  employ  without  actually 
solving  the  complete  mathematical  problem  and  comparing  the  results  to  measurements. 

In  this  report,  we  consider  the  simplest  approximate  form. 

To  obtain  the  approximate  form  we  shall  use,  we  first  evaluate  the  contribution  to 
the  buoyancy  force  due  to  mixing  with  V  •  D  given  by  Eq.  (22).  From  Eq.  (16)  the 
contribution  to  the  buoyancy  force  is 

d0*  [O  <  r  <  r(t),  0  <  t  <  tm 

0— 7—  (z0  -  z)9  for  J 

dz  I  0  <  r  <  r(tm ),  '  ^ 


t  >  tn 


and  zero  elsewhere. 
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Figure  4  shows  the  distribution  of  this  force.  The  net  downward  force  acting  on  the 
upper  half  of  the  region  in  Fig.  4  is 


dp* 


3  _3 


rot 


for  0  <  f  <  f_ 


dp* 

dz 


v  _  2  „  dP*  43  3. 

m)  '  3  ^  dz  A  °f" 


for  <  >  t 


m  • 


(24) 


By  symmetry,  a  force  equal  in  magnitude  acts  upward  on  the  lower  half  of  the  region. 


4- 

i  t 


A  * 


‘*4 


BOUNDARY  OF  MIXED  REGION  r(t) 


Fig.  4— Form  of  contribution  to  the  buoyancy  force 
due  to  mixing.  Arrows  indicate  relative  magnitude 
and  direction  of  force.  The  force  is  constant  along 
constant-;  lines  inside  the  region  but  changes  with 
vertical  distance  from  2  >  20;  outside  the  region  the 
force  is  zero. 


We  now  develop  an  alternate  form  of  V  •  D  that  produces  the  same  net  force.  First 
we  take  V  *  D  to  be  nonzero  over  a  square  region  of  dimension  2C;  a  square  region  is 
used  for  mathematical  convenience  in  later  calculations.  Within  this  square  we  take 

-V  •  D  -  (*0  -  z)  j-  for  0  <  f  <  fm,  (25) 

az  ‘m 

and  zero  otherwise.  This  leads  to  a  contribution  to  the  buoyancy  force  from  mixing  of 
the  form 


and 


a  dp*  1 
<*» 


for  0  <  f  <  tm 


-  zrt 


for  f  >  tm 


1 

I 

\ 

i 


(26) 
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inside  the  square  region.  The  net  buoyancy  force  due  to  mixing  acting  on  the  upper  half 
of  the  square  region  is 

0g  g3  7—  for  0  <  t  <  tm 

d*  lm 

for  t  >  tm.  (27) 

az 

Comparing  Eq.  (27)  with  Eq.  (24),  we  see  that  the  net  force  will  agree  at  each  instant  of 
time  if  we  take 

c  =  /f  r<'m)  =  °-874^^-  (28) 

We  note  that  for  t  >  tm  the  detailed  force  distribution  given  by  Eq.  (26)  has  precisely 
the  same  behavior  as  in  Eq.  (23),  although  the  force  now  acts  on  a  slightly  different  area. 
Equations  (25),  (17),  and  (28)  specify  the  forcing  function  we  shall  use  in  our  evaluation 
of  the  interned  waves  produced  by  the  growth  and  collapse  of  a  mixed  region. 


Solution  for  Fluid  with  Constant  N  and  Constant  Depth 

In  obtaining  a  mathematical  solution  to  the  problem  we  shall  first  solve  Eqs.  (10), 
(11),  and  (12)  in  terms  of  the  divergence  of  the  mean  turbulent  density  flux  V  •  D; 
afterward  the  solution  will  be  made  specific  for  the  mean  density  flux  described  by  Eq. 
(25).  This  procedure  reveals  the  necessity  of  employing  simple  analytical  forms  for  the 
divergence  of  the  mean  density  flux  if  the  numerical  problem  of  evaluating  the  disturb¬ 
ance  is  to  be  practicable. 

We  assume  at  the  outset  that  V  •  D  vanishes  everywhere  outside  a  square  region  of 
side  2fi  centered  on  the  2  axis  at  a  depth  z0\  furthermore,  V  •  D  is  taken  to  be  an  even 
function  of  x.  We  consider  the  fluid  to  have  a  constant  depth  D  with  constant  Brunt- 
Vaisala  frequency  N  given  by 


N2  =  J_  *£l 

Po  dz 

The  horizontal  extent  of  the  fluid  is  taken  as  infinite  (if  vertical  boundaries  are  present 
they  are  regarded  as  nonreflecting).  As  already  mentioned,  we  deal  with  situations  in 
which  the  mean  motion  is  two  dimensional,  dependent  on  x,  z,  and  t.  Figure  5  shows 
the  geometric  situation.  The  initial  conditions  are  that 


v  =  p=  0atf  =  0, 


and  the  boundary  conditions  are 


w  =  0  on  z  =  0,  D 
w(x,  z,  f)  remains  bounded  as  x~*  ±  00 ■ 


(29) 


(30) 
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Fig.  5— Location  of  region  in  which  mean  turbulent  density  flux 
0  is  nonzero.  The  density  p*(z)  is  the  initial,  linear  stratification. 


We  introduce  the  Fourier  transform  in  time  defined  by 

w(x,  z,  to)  =  I  w(x,z,  t)e~,utdcj, 


(31) 


with  similar  expressions  for  the  other  dependent  variables.  The  transformed  quantities 
will  be  analytic  functions  of  u>  in  some  region  Im(io)  <  cj0  (cj0  positive)  of  the  complex 
oj  plane.  From  the  inversion  theorem  we  get 


w(x 


1 

,z,t)  =  —  J  w(x,z,a})e,utdcj. 


.oo_,e 

'■-oa-ie 

with  e  chosen  so  that  the  path  of  integration  lies  below  the  singularities  of  w. 


(32) 


Applying  the  transform  to  Eqs.  (10)-(12)  and  using  the  initial  conditions  in  Eq.  (29), 
we  find  after  some  elimination  that 


w  + 


(33) 


in  which  subscripts  denote  differentiation,  and 

foo 

(-V  •  D)e-'wtdw. 

u 

The  boundary  conditions  for  w  follow  from  Eq.  (30)  and  can  be  written  as 

u>  -  0  on  z  =  0,D 
w  remains  bounded  as  x  -*  ±  00 . 


(34) 


(35) 


We  now  want  to  find  w  ( x ,  z,  u>). 
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Considering  w  (x,  z,  co)  as  a  function  of  z,  we  can  expand  it  in  a  Fourier  series 

OO 

U)(x,Z,U })  =  12  <t>n(z)Fn(x,  co). 


in  which 


*„(*,-  sin  (^f-) 


and  the  Fn  coefficients  depend  on  x  and  go.  The  motivation  for  using  the  orthogonal 
functions  sin  (ntrz/D)  rather  than  some  other  set  of  orthogonal  functions  follows  from 
the  fact  that  the  trigonometric  functions  are  eigenfunctions  of  the  differential  operator 
appearing  in  Gq.  (33). 

Substitution  of  the  expansion  given  by  Eq.  (36)  into  Eq.  (33)  and  use  of  Eq.  (37) 
and  the  orthogonal  property  of  produces  for  the  determination  of  Fn 

& F-  *  5,  f  k  (3 


in  which 


We  define 


a2  =  /^IZL\2 _ _ 

"  \D  )  N2  -  co2 

/"(JC’W)  =  N2  -  Cl)2  jf  ^fdZ' 


so  that  Eq.  (38)  can  be  written  in  the  form 

F  +  a2F  =  d2-  -  ■  (41) 

dx2  "  n  n  dx 2  ’  1  ' 

the  boundary  conditions  for  Fn  are  that  it  remain  bounded  as  x  -*  ±  °°. 

We  now  wish  to  find  the  solution  F„  (x,  co)  to  Eq.  (41).  The  dependence  of  Fn  on 
co  is  affected  by  an  which  also  depends  on  co.  From  Eq.  (39),  which  gives  o^(co)  we 
define  a  single- valued  root  as 


«n(w)  =  — 


'  D  y/N 2  -  co2  ’ 

in  which  we  use  the  principal  branch  for  y/N 2  -  co2  .  In  the  complex  co  plane,  then, 
branch  cuts  occur  for  Im  (co)=  0,  icoi  >  an(co)  is  not  defined  on  these  branch  cuts, 
but  it  possesses  values  through  out  the  remainder  of  the  plane  as  given  by  Eq.  (42). 
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Some  properties  of  an  needed  in  the  construction  of  the  solution  for  Fn  are  that  the 
imaginary  part  of  an  is  negative  for  1m  (co)  <  0  and  equal  to  zero  on  the  line  segment 
Im  (cj)  =  0,  |cu|  <  N. 

The  solution  to  Eq.  (41)  follows  from  the  variation-of-parameters  technique,  in 
which  we  set 


Fn(x)  =  an(x)eia-x  +  bn(x)e 

and  use  Eq.  (41)  with  its  boundary  conditions  to  find  an (.t )  and  bn(x).  This  process 
yields 


la. 


,ianX  J  g-'Wn 


(r?kfT7  +  e  <Xf 


eia"^ln  (17)4/17 


(43) 


It  is  easy  to  show  that  since  j'  (jc  ,  z,  a>)  is  an  even  function  of  x,  then  Fn  possesses  the 
same  property. 

We  now  have  the  general  form  of  the  solution  for  the  vertical  velocity  component  as 

(44) 


—  /  v  r°°  ‘  ^ 

w(x,z,t )  =  ^sinf^p);^-  f  Fn(x,cj)eiuldu, 

'  1  J  oo-  ,e 

with  Fn  given  by  Eq.  (43)  and  the  path  of  integration  located  below  the  singularities  of 
Fn.  Since  Fn  is  an  even  function  of  x,  we  conclude  that  the  vertical  velocity  compo¬ 
nent— in  fact  the  complete  disturbance— is  an  even  function  of  x,  as  we  would  expect 
from  the  physics  of  the  problem.  We  can  therefore  confine  our  attention  to  the  region 
x  >  0. 


Equation  (44)  gives  the  form  of  the  solution  when  the  only  assumptions  on  f  are 
that  it  be  an  even  function  of  x  vanishing  outside  some  finite  region.  Four  integrations 
occur  in  the  solution :  one  to  obtain  the  function  f ,  a  second  for  /„ ,  a  third  for  Fn , 
and  the  fourth  for  the  vertical  velocity  component  w  as  exhibited  explicitly  in  Eq.  (44). 
For  the  sake  of  convenience  in  performing  the  numerical  evaluation  of  the  disturbance, 
it  is  clearly  desirable  to  be  able  to  analytically  perform  as  many  of  these  integrations  as 
possible.  This  was  our  motivation  for  selecting  the  function  given  in  Eq.  (25)  in  place  of 
that  given  in  Eq.  (22). 

We  now  develop  the  specific  form  of  the  solution  when  the  mean  turbulent  density 
flux  satisfies  Eq.  (25).  From  Eq.  (34)  we  find  that 

f  =  0^-Uo  -  z)  —  7"<~'C0tm  *  ^  -  yWe  +  .v)  ,  (45) 

where  we  restrict  our  attention  to  the  region  x  >  0;  Eq.  (28)  specifies  the  quantity  8. 

This  result  yields  for  In  defined  in  Eq.  (40), 
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N2  1  -  exp  (~iutm  ) 

CO7 


N*~ 


//(«  -  x). 


(46) 


We  write  /„  as 


/„  =  CnH( 8  -  x)  (47) 

in  which  Cn  can  be  read  off  from  Eq.  (46).  Cn  is  independent  of  x,  and  Eq.  (47)  will 
allow  Fn  to  be  written  in  a  condensed  form. 

From  Eqs.  (43)  and  (47)  we  find  that 

Fn  =  Cn  {cos  (anx)e~ian^H(i  -  x)  -  ie'iQnX  sin  (a„C)//(x  -  «)}  ,  (48) 

with  a„  defined  by  Eq.  (42).  As  a  function  of  uj,  Fn  has  singularities  at  uj  =  ±  N,  and 
the  branch  cuts  for  an  leave  Fn  undefined  on  these  cuts.  However,  Fn  is  an  analytic 
function  of  to  below  the  real  axis.  In  Eq.  (44),  then,  e  can  be  taken  as  any  positive  num¬ 
ber.  Application  of  Cauchy’s  theorem  permits  us  to  transfer  the  integration  path  in  Eq. 

(44)  to  an  integration  along  the  real  axis.  In  carrying  out  this  process  we  must  take 
proper  account  of  the  singular  points  and  also  define  Fn  along  the  branch  ruts  through 
a  limiting  process. 

Substitution  of  Eq.  (48)  into  Eq.  (44)  and  use  of  the  procedure  just  outlined  yields 
the  specific  form  of  the  vertical  velocity  component  w  for  the  model  under  consideration. 
Different  expressions  exist  for  0  <  x  <  8  and  x  >  H .  For  both  regions  we  find  that  we  can  write 

u>(x,z,t)  =  w_{x,z,t)  +  u>+(x,z,  t),  (49) 

in  which  w_  involves  an  integration  over  frequencies  0  ^  o>  <  N  and  ie+  an  integration 
over  frequencies  cj  >  N.  These  two  functions  have  the  following  explicit  forms.  For 
0  <  x  <  C, 


u>.(*,z,  f) 


[sin  (cut  -  crng)  -  sin  (cut  -  u )tm  -  on6)l 

w*m  J 


(50) 


in  which 
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(51) 


and 


/  4  1  .  ( nitz\  ( nirz0  \  \(  nil  g  \  (m rS  \  (niri\ 

^sm(-D-  cos{-D^)  [  D-jcos  (t)  "  81n  \~F / 

fl-1  *- 


r  n2 

l  exp  ,WJ 


7,8)  cosh  (0„x)  - 


sin  ut  -  sin  a )(t  - 


ut. 


■^1  d«, 

J  (! 


in  which 


mr  u) 

~  D  -  N2  ‘ 


(52) 


(53) 


For  the  region  x  >  6,  we  have 


w . 


(  nnz  > 

\  cos  /  22*ii  \ 

(nit 8  \ 

(—) 

/ 

\  D  i 

IC08(dJ 

ssr) cos 

/ 

-  sin  1 

rv  n2 

X  1  N2  -  M«  ,ln 


“0 
X  dcj, 

in  which  an  is  as  given  in  Eq.  (51),  and 


in  (a„8)  £ 


cos  (u)t  -  u)tm  -  a„x)  -  cos(a;t  -  anx) 


(54) 


£  »"  (if)  l2?0)  [(^) CM  (if)  - sin  (tf) 

n-1  I-  J 

r00  ^2  sin  u(t  -  f  _  )  -  sin  u>fl  , 

X  A/2  _  w2  exP  (-^«*)  8Ulh  <M>  [ - ^ - J  <*<*>,  (5B> 

in  which  /J„  is  as  in  Eq.  (53). 

The  above  expressions  provide  the  vertical  velocity  component  throughout  the  fluid 
as  a  function  of  time.  We  can  easily  determine  the  remaining  quantities  describing  the 
disturbance  from  Eqs.  (10M13).  In  particular  the  displacement  t}  defined  by  Eq.  (13) 
follows  from  an  integration  over  time  of  the  expressions  in  Eqs.  (50),  (52),  (54),  and 
(55);  since  the  results  can  be  written  immediately  from  the  cited  equations,  we  shall  not 
write  out  the  explicit  forms.  In  our  later  comparison  to  experimental  results  we  use  the 
temperature  change  AT  at  a  point  outside  the  mixed  region,  which  can  be  expressed  in 
terms  of  the  displacement  by  use  of  Eq.  (14). 
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We  note  that  all  quantities  describing  the  disturbance  can  be  written  as  the  sum  of 
two  terms;  the  first  one  contains  an  integration  over  frequencies  0  <  co  <  N,  while  the 
second  has  an  integration  over  frequencies  w  >  N.  The  term  involving  frequencies  less 
than  the  Brunt-Vaisala  frequency  N  represents  a  propagating  system  of  internal  waves. 
The  term  with  frequencies  greater  than  N  does  not  represent  propagating  waves;  rather,  it 
is  a  forced  disturbance  maintained  by  the  mean  turbulent  density  flux.  This  forced  dis¬ 
turbance  decays  exponentially  with  horizontal  distance  from  the  boundary  of  the  mixed 
region  and  decays  at  sufficiently  large  times  as  well. 


COMPARISON  WITH  EXPERIMENTAL  RESULTS 

Before  proceeding  with  our  comparison,  we  make  some  comments  on  the  degree  of 
mixing  0.  Direct  measurements  of  0  do  not  exist.  We  used  0  =  0.5  in  our  calculations 
and  cite  some  indirect  experimental  evidence  supporting  this  choice.  In  Fig.  2a,  if  we 
assume  that  at  the  time  of  maximum  vertical  expansion  a  constant  degree  of  mixing  0 
exists  in  the  mixed  region,  then  we  have 

.  final  vertical  dimension  __ 

u  —  —  0.55. 

maximum  vertical  dimension 

Similarly,  Schooley  and  Stewart  (1)  present  data  that  imply  a  value  0  =  0.60  at  the  time 
of  maximum  vertical  extent;  in  this  case  the  mixing  device  is  a  self-propelled  body  that 
moves  through  the  fluid.  These  approximate  values  for  0  apply  to  the  mixed  region  at 
the  time  of  its  maximum  vertical  extent.  However,  we  take  these  values  as  indicating 
that  various  mixing  devices  produce  approximately  the  same  degree  of  mixing  and  that  a 
constant  value  of  0  =  0.5  is  acceptable  as  a  first  approximation. 


Comparison  for  Large  Horizontal  Distances  from  the  Mixed  Region 

Schooley  and  Hughes  (9)  present  measurements  of  the  temperature  change  produced 
at  points  near  the  top  surface  of  a  stratified  fluid  as  a  consequence  of  localized  mixing  in 
the  interior  of  the  fluid.  Figure  6  shows  their  measured  data  for  two  stratifications, 
which  can  be  approximated  by  constant  values  of  the  Brunt-Vaisala  frequency  N.  At 
least  two  sets  of  measured  data  were  obtained  for  each  stratification.  Pertinent  values  of 
the  parameters  are  listed  below. 


N  (s’1 ) 

dT*ldz,  (°C/cm) 

z  (cm) 

1.19 

11.1 

0.250 

2.0 

0.68 

3.81 

0.300 

2.0 

D  -  7.3  cm, 

Zq  =4.5  cm. 

r0  =  0.65  cm 

The  quantity  z  gives  the  depth  at  which  measurements  were  made;  for  each  stratification 
temperature  as  a  function  of  time  was  measured  at  the  three  horizontal  distances 
x  *  5.72,  8.26,  and  12.7  cm. 

The  solid  curves  in  Fig.  6  are  from  the  theory  of  Schooley  and  Hughes  (9).  Their 
theory  assumes  that  wave  generation  does  not  begin  until  the  mixed  region  has  attained 
its  maximum  vertical  dimension;  furthermore,  the  mixing  is  assumed  complete  so  that 


* 


►— I04°C  I - 1 0  2°C 

N  =  I  l»Mt"  N  =0  68  we'1 


Fig.  6— Measured  and  theoretical  values  of  temperature  changes  produced  by  growth  -.nd 
collapse  of  a  mixed  region. 

. measured  data  of  Schooley  and  Hughes  (9), - theory  of  Schooley  and  Hughes 

(9),  ©  according  to  theory  developed  in  this  report  Values  at  the  lop  of  each  curve  give 
the  horizontal  distance  x  in  centimeters  at  which  measurements  and  calculations  were  made, 
(based  on  Schooley  and  Hughes  (9)) 


0=1.  Thus,  their  mathematical  model  represents  wave  generation  as  the  result  of 
collapse  of  an  initial  density  perturbation  in  a  stratified  fluid  initially  at  rest.  They  take 
their  time  origin  at  the  moment  of  collapse.  In  comparing  their  theoretical  values  with 
the  measured  values  they  shifted  their  curves  until  the  best  fit  was  obtained;  this  was 
considered  permissible  since  the  actual  moment  of  collapse  was  unknown.  Figure  6 
shows  that  their  theory  yields  fair  agreement  with  amplitudes  but  appears  to  indicate 
more  oscillations  than  are  observed. 

Our  theory  assumes  that  the  mixed  region  expands  from  a  point  according  to  the 
rule  given  by  Eq.  (19).  A  value  is  needed  for  the  quantity  A  which  appears  in  Eqs.  (19) 
and  (28).  An  approximate  value  of  A  =  2.06  s_,/3  was  obtained  from  the  measured  radii 
of  the  mixed  region  at  the  time  of  collapse  by  setting  the  measured  radius  equal  to 
Ar0(/m  ),/3  f°r  each  stratification  and  then  averaging  the  resulting  values  of  A.  Values  of 
the  measured  radii  are  given  by  Schooley  and  Hughes  (9). 

The  circled  points  in  Fig.  6  were  calculated  with  the  theory  developed  in  this  report. 
In  evaluating  the  modal  sum  we  replaced  the  sum  over  all  the  modes  by  a  partial  sum 
containing  a  sufficient  number  of  modes  to  provide  a  value  for  the  displacement  r? 
stable  to  ±  0.001  cm.  For  both  stratifications  the  number  of  modes  required  increased 
with  time  at  each  location.  A  minimum  number  of  four  modes  was  calculated  for  each 
point.  The  maximum  number  of  modes  required  decreased  with  horizontal  distance  from 
the  mixed  region.  Thus  for  N  =  1.19  s*1  the  maximum  number  of  modes  was  30  at  5.72 
cm,  21  at  8.26  cm,  and  14  at  12.7  cm.  Approximately  2/3  as  many  modes  were  required 
at  the  same  locations  for  N  =  0.68  s“l . 

On  the  whole,  our  theoretical  values  show  a  fair  agreement  with  the  measured  ampli¬ 
tudes  and  phases.  We  should  point  out  that  Schooley  and  Hughes  made  their  measure¬ 
ments  very  close  to  a  rigid  boundary  and  that  they  included  in  their  theory  a  correction 
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for  the  viscous  layer  at  the  boundary.  It  seems  that  if  our  theory  were  modified  to 
account  for  the  boundary  layer  effect,  we  would  improve  our  agreement  with  the  meas¬ 
ured  data  because  our  calculated  amplitudes  would  be  reduced  at  later  times. 

Although  our  theoretical  values  also  show  some  agreement  with  the  theory  of 
Schooley  and  Hughes,  the  basic  premises  of  the  two  theories  are  quite  different.  Our 
model  describes  the  disturbance  as  originating  at  the  start  of  mixing  instead  of  at  the 
time  of  maximum  vertical  expansion;  thus  there  is  no  uncertainty  about  the  meaning  of 
our  time  origin,  and  the  model  accounts  for  the  entire  growth  phase  of  the  mixed  region. 
Furthermore,  we  use  a  constant  degree  of  mixing,  0  =  1/2,  which  seems  partially  sup¬ 
ported  by  experimental  evidence,  while  the  calculations  of  Schooley  and  Hughes  assume 
that  0-1. 


Comparison  Near  the  Mixed  Region 

We  now  consider  the  experimental  results  of  Williams  (6).  He  made  temperature 
measurements  at  a  point  just  outside  the  mixed  region  at  its  maximum  vertical  extent. 
These  measurements  describe  the  near-field  behavior  of  the  disturbance  caused  by  mixing 
in  contrast  to  the  far-field  behavior  measured  by  Schooley  and  Hughes  (9).  Moreover, 
Williams  studied  the  effect  of  excitation  time  te  (the  time  interval  during  which  the  mixer 
agitates  the  fluid)  on  the  disturbance. 

He  used  an  axisymmetric  mixer  consisting  of  a  rod  wound  with  wire  having  a  total 
effective  radius  r0  =  0.95  cm.  The  mixer  was  placed  at  a  depth  z0  =  9.5  cm  in  fluid  of 
total  depth  D  -  25.4  cm.  Measurements  were  made  at  x  =  1.35  cm,  z  =  7.5  cm.  The 
stratification  of  the  fluid  onsisted  of  two  approximately  homogeneous  layers  at  the  top 
and  bottom  with  an  intermediate  layer  having  an  approximately  linear  density  variation. 
Excitation  time  te  varied  from  1  s  to  32  s. 

Before  proceeding  with  a  comparison  of  our  theory  with  the  measurements  we  men¬ 
tion  several  reasons  why  our  theory  is  not  entirely  applicable  to  the  experiment. 

1.  Our  theory  assumes  that  a  small  disturbance  exists  in  the  fluid,  and  this 
assumption  may  fail  to  hold  near  the  mixer. 

2.  We  used  a  simplified  expression  for  the  divergence  of  the  turbulent  density  flux 
V  •  D  in  order  to  evaluate  several  integrals  that  appear  in  the  theory.  This  simplified  ex¬ 
pression  may  produce  a  poor  approximation  to  the  near-field  behavior  but  a  good  approx¬ 
imation  to  the  far-field  behavior. 

3.  The  presence  of  the  solid  mixer  requires  the  fluid  to  flow  around  it,  and  this 
was  not  accounted  for  in  our  boundary  conditions. 

4.  Stratification  of  the  fluid  used  in  the  experiment  does  not  follow  the  linear 
behavior  (constant  N)  assumed  by  the  theory. 

5.  The  measurements  apply  to  a  particular  realization  of  the  flow.  On  the  other 
hand,  the  theory  deals  with  the  ensemble  average  disturbance.  It  is  not  clear  to  what 
extent  particular  realizations  deviate  from  the  ensemble  average  outside  the  mixed  region. 
This  matter  will  eventually  have  to  be  settled  by  experiment.  Several  of  these  effects 
seem  quite  important  for  our  comparison  because  the  near-field  behavior  intimately 
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depends  on  the  detailed  flow  processes  occurring  in  the  mixed  region.  Despite  these 
qualifications  we  are  motivated  to  make  a  comparison  by  the  belief  that  some  of  the 
gross  features  contained  in  the  measured  data  should  also  appear  in  the  theoretical  results. 


In  order  to  apply  our  theory  to  the  experiment  we  must  determine  a  value  for  the 
parameter  A  entering  Eq.  (18).  This  parameter  is  related  to  the  size  of  the  mixed  region 
Rather  than  assume  a  virtual  origin  for  the  mixed  region,  we  shall  use  Eq.  (18)  in  which 
t  *  0  denotes  the  start  of  mixing.  Thus  at  t  =  0  we  require  that  the  boundary  of  the 
mixed  region  coincide  with  the  effective  radius  of  the  mixer;  this  gives 


r(0)  =  r0  =  Aroyt^  , 


which  establishes  a  relation  between  A  and  t0 .  The  solid  curves  in  Fig.  7  show  the  meas¬ 
ured  temperature  oscillation  at  a  distance  2.42  cm  from  the  center  of  the  mixer.  Turbu¬ 
lent  fluctuations  do  not  appear.  However,  for  the  greatest  excitation  time  of  32  s  there 
are  high-frequency  oscillations  which  we  take  as  an  indication  that  the  boundary  of  the 
mixed  region  is  in  the  vicinity  of  the  measured  point.  We  find  from  Eq.  (18)  that 
A  as  0.7  s-1'*  by  requiring  the  maximum  radius  of  the  mixed  region  r(tm )  to  equal  the 
distance  of  the  measurement  point  from  the  center  of  the  mixer.  This  is  the  approximate 
value  we  employed  in  our  calculations. 


i 
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Fig.  7-Temper*  ure  oscillation  at  a  point  for  various  excitation  times.  Conditions  are  D  *  25.4  cm, 
*0  ■  9.5  cm,  rj)  ■  0.95  cm,  x  ■  1.35  cm,  t  *  7.5  cm.  The  Brunt- Vaisala  period  increases  from  toy  •  26.6s 
in  (a)  to  31  s  in  (f) 


measurements  of  Williams  (6), . theory 
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We  account  for  the  presence  of  the  mixer  in  the  following  indirect  fashion.  We  con¬ 
sider  the  mixed  fluid  to  lie  between  the  outside  radius  given  by  Eq.  (18)  and  the  surface 
of  the  mixer,  which  has  an  effective  radius  r0 .  In  this  region  the  turbulent  density  flux 
is  assumed  to  satisfy  Eq.  (22)  in  which  t  =  0  now  designates  the  start  of  mixing.  Using 
the  procedure  following  Eq.  (22)  we  again  are  led  to  the  approximation  for  the  turbulent 
density  flux  given  by  Eqs.  (25)  and  (28).  This  scheme  partially  accounts  for  the  fact  that 
mixed  fluid  does  not  fill  the  region  occupied  by  the  mixer. 


Figure  7  shows  a  comparison  between  measured  temperature  changes  (solid  curves) 
and  temperature  changes  according  to  our  theory  for  100  s.  (We  note  that  for  times  in 
excess  of  120  s  there  may  be  interference  in  the  measured  data  due  to  wave  reflections 
from  the  walls  of  the  tank.)  While  the  comparison  shows  no  agreement  in  detail,  some 
general  properties  agree.  Our  theory  gives  an  increase  in  temperature  as  soon  as  mixing 
begins,  indicative  of  a  downward  motion  of  the  fluid.  This  agrees  qualitatively  with  the 
measurements,  but  our  calculated  temperature  does  not  demonstrate  the  rapid  initial 
increase  found  in  the  measurements. 


A  further  point  of  qualitative  agreement  arises  when  we  compare  the  behavior  of 
average  properties  of  the  temperature  fluctuation.  The  properties  we  compare  are  the 
average  temperature  fluctuation  A H  and  the  root-mean -square  (rms)  temperature  deviation 
o  from  the  average.  For  the  measured  data  these  averages  were  calculated  for  the  time 
range  0  <  t  <  120  s  by  using  values  of  the  temperature  fluctuation  at  4-s  intervals.  For 
the  theoretical  curves  the  averages  were  calculated  for  the  time  interval  0  <  t  <  100  s  at 
5-s  intervals.  Thus,  AW  and  a  are  defined  by 


and 


o  =  -  AH)2  , 

in  which  A Tn  is  the  temperature  change  for  the  interval  n.  Figure  8  shows  how  these 
two  quantities  vary  with  excitation  time  according  to  the  measured  data  and  according 
to  theory.  All  the  curves  demonstrate  a  trend  to  increase  with  excitation  time,  which 
roughly  indicates  that  the  disturbance  grows  in  magnitude  with  an  increase  of  excitation 
time. 


Figure  9  shows  the  displacement  as  a  function  of  time  at  various  locations  in  the 
fluid  for  an  excitation  time  te  =  5  s.  These  curves  were  calculated  using  the  theory  to 
demonstrate  how  the  features  of  the  signal  change  with  distance  from  the  mixer  and  with 
depth  in  the  fluid.  An  interesting  feature  shown  by  the  curves  is  that  the  temperature 
fluctuation  has  a  more  complex  form  in  the  interior  of  the  fluid  (Fig.  7a  and  7b)  than  it 
has  near  the  surface  (Fig.  7c). 

The  comparison  seems  to  indicate  that  our  theory  is  inadequate  near  the  mixed 
region.  The  major  reason  for  this  inadequacy  may  be  that  our  theory  does  not  lead  to  a 
disturbance  increasing  in  amplitude  at  a  sufficient  rate  during  the  initial  phases  of  the 
motion,  as  demonstrated  in  Fig.  7.  We  base  this  remark  on  the  expectation  that,  if  the 
first  maximum  temperature  deviation  were  larger,  then  the  amplitudes  of  later  tempera¬ 
ture  extremes  would  be  correspondingly  increased.  The  fact  that  our  theory  does  not 
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Fig.  8— Average  characteristics  of  the  temperature  oscillations  as  a 
function  of  excitation  time 

- measurements  of  Williams  (6), . theory 

produce  a  sufficiently  large  initial  disturbance  must  of  course  be  associated  with  our 
forcing  function  V  •  D. 

We  recall  that  we  replaced  the  expression  for  V  •  D  in  Eq.  (22)  by  the  expression 
in  Eq.  (25)  on  the  grounds  of  mathematical  convenience.  Actually,  V  •  D  given  by  Eq. 
(25)  corresponds  to  smearing  out  V  •  D  in  Eq.  (22)  over  a  prescribed  area.  The  buoy¬ 
ancy  force  due  to  mixing  is  also  spread  out  over  the  prescribed  area  but  such  as  to  pre¬ 
serve  the  net  force  acting  on  the  fluid.  What  this  means  is  that  at  most  points  in  the 
region  of  mixing  we  weakened  the  buoyancy  force  due  to  mixing.  It  seems  reasonable  to 
expect  that  the  disturbance  at  a  point  near  the  mixed  region  will  depend  more  on  the 
actual  distribution  of  the  buoyancy  force  due  to  mixing  than  on  its  intergral  properties. 
Thus,  it  may  be  that  in  smearing  out  the  buoyancy  force  we  weakened  the  disturbance 
in  the  near  field. 

We  give  some  suggestions  later  as  to  how  the  present  theory  may  be  improved  by  the 
use  of  an  expression  for  V  •  D  which  should  approximate  V  •  D  in  Eq.  (22)  in  such  a  way 
as  to  avoid  the  smearing  out  process,  and,  consequently,  the  weakening  of  the  near-field 
disturbance. 


CONCLUSION 

We  developed  in  this  report  an  approximate  mathematical  model  to  describe  the  in¬ 
ternal  wave  disturbance  caused  by  the  growth  and  collapse  of  a  localized,  two-dimensional 
region  of  turbulence  in  a  linearly  stratified  fluid.  Our  theory  differs  from  earlier  work  to 
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RELATIVE  LOCATIONS  AT  WHICH 
CALCULATIONS  MADE  (VERTICAL 
SCALE  •  1/2  HORIZONTAL  SCALE) 


Fig.  9— Displacement  as  a  function  of  time  at  various  positions  in  the  fluid  (according  to  theory). 
The  conditions  are  zq  -  9.5  cm,  te  -  5  sec,  dT*ldt  m  ,3510°C/cm,  tgy  »  27.7  sec;  in  (a),  (b)and 
(c)  the  curves  from  bottom  to  top  are  for  lateral  distances  x  ■  20,  40,  and  80  cm  respectively. 


large  distances  from  the  mixed  region;  however,  comparison  with  data  describing  the  dis¬ 
turbance  near  the  mixed  region  indicated  that  our  model  does  not  provide  an  entirely 
faithful  representation  close  to  the  mixed  region.  We  present  further  details  of  the  com¬ 
parison  after  we  summarize  the  essential  features  of  the  model. 
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Our  model  applies  to  the  mean  disturbance,  determined  as  an  ensemble  average  of 
the  motion.  We  linearized  the  Navier-Stokes  equations  in  the  form  applicable  to  the 
mean  disturbance  to  obtain  equations  amenable  to  mathematical  analysis  whose  applica¬ 
bility  to  the  physical  process  could  then  be  studied.  The  forcing  function  for  the  dis¬ 
turbance  in  the  linearized  system  of  differential  equations  is  the  divergence  of  the  mean 
density  flux  V  •  D  associated  with  the  turbulent  motion.  This  density  flux  describes  how 
the  mean  density  in  the  turbulent  region  changes  as  a  consequence  of  mixing.  Specifica¬ 
tion  of  V  *  D  in  functional  form  constituted  a  basic  problem  of  the  model. 

Intuitive  and  empirical  considerations  provided  an  analytical  expression  for  V  •  D  for 
mixed  regions  which  can  be  described  by  a  constant  degree  of  mixing  0.  Further  assump¬ 
tions  that  entered  the  argument  establishing  the  expression  for  V  •  D  were 

1.  Mixing  in  the  fluid  occurs  during  a  finite  time  interval  of  length  tm  -  0.2 tBV  + 
te  in  which  tB  v  denotes  the  Brunt-Vaisala  period  and  te  the  time  interval  during  which 
the  mixer  device  agitates  the  fluid. 

2.  The  mixed  region  can  be  regarded  as  circular  during  the  time  interval  0  <  f  < 
tm  (t  =  0  represents  the  virtual  origin  of  the  mixed  region)  with  a  radius  increasing  ac¬ 
cording  to  a  t1/3  law. 

3.  p  =  1/2  represents  a  reasonable  value  for  the  degree  of  mixing  produced  by 
mixer  devices. 

We  cited  limited  experimental  support  for  some  of  these  assumptions. 

It  was  pointed  out  that  the  above  model  leads  to  a  rather  complicated  numerical 
problem  due  to  the  presence  of  several  integrals  in  the  solution.  For  the  sake  of  simpli¬ 
fying  the  numerical  problem  we  replaced  the  analytical  expression  for  the  flux  term 
V  •  D  by  an  alternate  expression  which  allowed  several  of  the  intergrals  to  be  evaluated 
explicitly.  Our  alternate  expression  for  V  •  D  is  equivalent  to  the  original  expression  in 
the  following  sense.  At  any  instant  of  time  the  two  expressions  lead  to  the  same  net 
buoyancy  force  associated  with  mixing  for  the  upper  (or  lower)  half  of  the  mixed  region. 

Using  this  alternate  expression  for  the  forcing  term,  we  obtained  a  solution  for  the 
vertical  velocity  component  in  the  case  of  a  linearly  stratified  fluid  of  constant  depth. 

The  other  quantities  describing  the  disturbance  can  be  immed  iately  evaluated  in  terms  of 
this  velocity  component.  In  particular,  the  vertical  displacement  of  a  particle  following 
the  mean  motion  is  the  time  integral  of  the  vertical  velocity  component.  The  importance 
of  the  vertical  displacement  stems  from  the  fact  that  outside  the  region  of  mixing  it  is 
directly  proportional  to  the  temperature  change  produced  by  the  disturbance,  which 
often  is  the  quantity  measured  in  experimental  investigations  of  internal  waves. 

We  compared  our  theoretical  solution  with  results  of  two  experimental  investigations. 
The  first  case  pertained  to  the  nature  of  the  disturbance  at  large  distances  from  the  mixed 
region.  Figure  6  shows  the  comparison  of  our  theory  with  measured  data  from  Schooley 
and  Hughes  (9)  for  two  values  of  the  Brunt-Vaisala  frequency,  and  with  the  results  of 
their  theory.  Our  theory  demonstrates  a  fair  agreement  with  the  measured  data.  We  note 
that  the  measurements  were  made  near  a  rigid  boundary  so  .hat  viscous  boundary-layer 
effects  were  probably  of  some  importance.  If  corrections  were  included  in  our  theory  to 
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account  for  the  viscous  layer  at  the  boundary  the  agreement  with  the  data  would  proba¬ 
bly  be  improved  because  the  disturbance  would  tend  to  be  diminished,  especially  at  later 
times. 

The  comparison  of  our  calculations  with  calculations  given  by  Schooley  and  Hughes 
(9),  shown  in  Fig.  6,  indicated  a  rather  close  agreement  during  the  early  stages  of  the 
motion.  However,  at  later  times,  our  theory  appears  to  yield  fewer  oscillations— in  agree¬ 
ment  with  the  measured  data — than  does  the  theory  of  Schooley  and  Hughes,  even  though 
our  theory  employs  premises  quite  different  from  those  used  by  Schooley  and  Hughes. 
For  example,  we  account  for  the  initial  growth  stage  of  the  mixed  region,  whereas 
Schooley  and  Hughes  assume  that  the  initial  growth  stage  can  be  ignored;  we  assume  a 
degree  of  mixing  0  =  1/2,  whereas  Schooley  and  Hughes  take  the  degree  of  mixing  0  =  1; 
in  our  theory,  internal  wave  generation  begins  at  the  same  time  as  does  mixing,  whereas 
Schooley  and  Hughes  assume  that  wave  generation  does  not  begin  until  the  mixed  region 
has  attained  its  maximum  vertical  extent  and  collapse  has  begun. 

The  second  case  pertained  to  measurements  made  at  a  point  very  close  to  the  mixed 
region  (Fig.  7).  While  agreement  in  detail  is  poor,  we  must  qualify  this  comparison  with 
the  remark  t>at  our  theory  is  not  directly  applicable  to  near-field  measurements  for 
specific  reasons  given  in  the  text. 

Thus  Fig.  7  shows  that  the  amplitudes  and  periods  are  significantly  different.  Since 
the  calculated  values  cannot  really  be  expected  to  agree  in  detail  with  near-field  measure¬ 
ments,  a  more  valid  comparison  can  be  made  on  the  basis  of  observed  trends  in  the  data. 
In  Fig.  8,  we  show  how  the  average  temperature  change  taken  over  a  finite  time  interval 
and  the  rms  temperature  change  from  its  average  over  the  same  time  interval  vary  with 
excitation  time.  The  trend  demonstrated  by  the  measured  data,  is  also  exhibited  by  the 
theoretical  values,  although  the  rates  of  increase  do  not  match  very  well. 

We  concluded  from  this  comparison  that  our  model  is  deficient  in  its  description  of 
the  disturbance  near  the  mixed  region.  A  partial  cause  of  this  deficiency  must  stem  from 
our  use  of  a  modified  expression  for  the  flux  term  V  •  D.  The  modified  expression  was 
used  in  order  to  simplify  some  of  the  mathematical  expressions  appearing  in  the  theory. 
However,  it  does  not  preserve  all  of  the  characteristics  of  the  original  expression  we 
found  for  V  •  D.  Since  the  disturbance  near  the  mixed  region  probably  depends  on  the 
detailed  nature  of  the  forcing  term,  we  must  expect  that  our  modified  expression  for 
V  •  D  does  not  give  a  good  description  of  the  disturbance  near  the  mixed  region. 

Some  remarks  follow  on  how  this  deficiency  might  be  corrected,  but  in  such  a 
manner  as  to  retain  mathematical  simplicity  of  the  theory. 


SUGGESTIONS  FOR  FURTHER  WORK 

We  pointed  out  that  our  approximation  to  the  forcing  term  V  •  D  given  in  Eq.  (22) 
by  the  expression  in  Eq.  (25)  is  probably  responsible  for  weakening  the  disturbance  in 
the  near  field.  We  now  outline  a  procedure  for  correcting  this  deficiency. 

We  retain  the  basic  premises  which  led  us  to  Eq.  (22).  The  problem  is  to  devise  an 
approximate  form  for  V  •  D  which  retains  the  essential  features  of  the  expression  in  Eq. 
(22)  and  at  the  same  time  allows  as  many  as  possible  of  the  integrals  appearing  in  the 
theory  to  be  evaluated  analytically. 


26 


J.  M.  BERGIN 


Since  V  *  D  vanishes  for  times  t  >  tm  we  only  have  to  consider  the  mixed  region 
for  the  interval  0  <  t  <  tm .  In  this  time  interval  we  suggest  that  (a)  the  boundaryt 
r(t)  =  Ar0j/7  of  the  mixed  region  be  replaced  by  a  square  region  of  the  same  area,  as  in 
Fig.  10a,  (b)  the  variation  of  the  square  boundary  with  time  then  be  approximated  by  a 
series  of  linear  laws,  as  in  Fig.  10b,  and  (c)  the  divergence  of  the  mean  turbulent  density 
flux  be  approximated  by 


-V  ■  D  =  (56) 

dz 

in  which  y  =  z0  -  z,  and  tb(x,  y )  designates  the  time  at  which  the  boundary  reaches  the 
point  (jc,  y).  This  expression  for  V  •  D  leads  to  a  density  change  caused  by  mixing  of 
the  form 


dp* 

Sp  =  p-^-izo  -  z) 


inside  the  square  region.  Thus  we  would  have  a  constant  degree  of  mixing  inside  the 
square  region. 


Fig.  10— (a)  Square  of  side  28(f)  used  to  replace  circular  region,  (b) - 8(f)  proportional  to  f1/3 

- series  of  linear  growth  laws  for  8(f)  to  approximate  the  solid  curve. 

(f  *  0  is  the  virtual  origin  of  the  mixed  region) 


In  Eq.  (56)  the  quantity  tb{x,  y)  depends  linearly  on  x  and  y  because  a  linear 
growth  rate  is  assumed  for  the  boundary.  This  property  of  tb(x,  y)  should  enable  the 
integrals  appearing  in  the  theory,  apart  from  the  final  integration  over  frequency,  to  be 
evaluated  in  analytic  terms  since  the  integrands  will  be  of  trigonometric  type. 


t  the  procedure  outlined  can  be  applied  to  any  growth  law  for  the  radius  of  the  mixed  region. 
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